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SUMMARY 
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This  report  presents  the  results  of  the  analysis  of  exceedance 
curves  by  aircraft  tail  number  for  the  F-4  and  A-37B  aircraft. 

The  aim  of  this  report  is  to  provide  a  statistical  bound  for 
the  different  exceedance  curves  that  are  obtained  from  various 
aircraft  tail  numbers  on  a  given  aircraft  type.  The  F-4  and  A-37B 
aircraft  were  selected  only  because  load  factor  (nz)  exceedance  data 
by  tail  number  were  available  in  table  formats  in  References  1  and  2. 

Results  of  this  report  indicate  that  reasonable  bounds  on 
exceedance  curves  may  be  obtained  by  tolerance  limits  which  assume 
that  for  a  fixed  load  factor  (nz),  the  log  of  the  exceedances  per 
thousand  flight  hours  has  a  normal  distribution.  Also,  tolerance 
limits  could  be  obtained  using  nonparametric  techniques  when  the 
quantity  of  tall  niimbers  is  sufficiently  large. 

Finally,  curves  were  fitted  to  the  composite  data.  The  designer 
may  extrapolate  these  curves  to  determine  design  limit  load  factor 
and  change  the  shape  of  the  curves  by  varying  values  of  coefficients. 
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DEFINITION  OF  SYMBOLS 


nz 

MED  (nz) 

Ei(nz) 

LCI(nz) 

UCI(nz) 

C 

MAX(nz) 

MIN(nz) 

P 

Y 

M(nz) 
s  (nz) 

NorUL(nz) 

No rLL(nz) 

M3=log^Q(C(nz)) 


vertical  load  factor  in  g's. 

the  median  value  of  exceedances  per  1,000  hours  for 
a  fixed  value  of  nz. 

the  number  of  exceedances  per  1,000  hours  for  a  fixed 
value  of  nz  for  the  ith  aircraft  tail  number. 

the  lower  bound  of  the  90%  confidence  interval  for 
MED(nz). 

the  upper  bound  of  the  90%  confidence  interval  for 
MED(nz). 

the  composite  exceedance  curve. 

the  maximum  number  of  exceedances  per  1,000  hours  for 
a  fixed  value  of  nz. 

the  minimum  number  of  exceedances  per  1,000  hours  for 
a  fixed  value  of  nz. 

the  percent  of  tall  numbers  that  lie  between  2  given 
limits . 

confidence,  i.e.  the  probability  that  a  given  statement 
is  true. 

mean  of  exceedances  per  1,000  hours  for  fixed  nz. 

standard  deviation  of  exceedances  pter  1,000  hours  for 
fixed  nz. 

upper  bound  of  90%  two-sided  tolerance  limit  with  90% 
confidence  assuming  normal  distribution. 

lower  bound  of  90%  two-sided  tolerance  limit  with  90% 
confidence  assuming  normal  distribution. 

used  as  the  mean  for  log  normal  assumptions. 


viii 


s' (nz) 

LNUL 

LULL 

LgNorUL 

LgNorLL 

F 

b(nz) 


DEFINITION  OF  SYMBOLS  (cont'd) 


standard  deviation  of  the  logs  of  the  Ei(nz)  for  a 
fixed  nz. 

Upper  bound  on  log  scale  of  90%  two-sided  tolerance 
limits  with  90%  confidence  of  exceedances  per  1,000 
hours . 

lower  bound  on  log  scale  of  90%  two-sided  tolerance 
limits  with  90%  confidence  of  exceedances  per  1,000 
hours . 

10 

(LNLL) 

curve  fitted  to  nz  composite  exceedance  per  4,000  hours, 
best  fit  polynomial  that  log  F  =  b(nz). 

Weighted  mean. 


SECTION  I 
INTRODUCTION 


As  a  result  of  a  series  of  meetings  with  Mr.  W.  J.  Crichlow, 
ASD/ENF,  it  appears  that  some  changes  are  necessary  in  the  pre¬ 
sentation  of  flight  loads  data.  This  is  expected  to  Improve  the 
development  of  more  definitive  structural  design  criteria.  Mr. 
Crichlow,  who  has  prime  responsibility  for  MIL-A-8866B  revision, 
explained  the  need  for  re-examining  the  statistical  variations  in 

vertical  load  factor  (n  )  exceedance  curves  and  this  distribution's 

z 

Impact  on  more  precise  aircraft  design,  and  fatigue,  and  fracture 

analyses.  Accordingly,  AFFDL/FBE  initiated  a  program  to  analyze 

flight  loads  data,  categorized  by  mission  segment,  type,  and  tail 

number  as  available,  to  determine  the  load's  distribution  mean 

exceedance  curve  and  the  statistical  spread  about  the  mean. 

Load  factor  (n  )  data  from  A-37B  (Reference  1)  and  F-4 
z 

(Reference  2)  aircraft  were  segregated  by  tail  number,  base,  and 
mission  type  as  available  in  the  reports.  The  data  from  these 
subdivisions  were  analyzed  and  the  results  plotted  in  the  form  of 
exceedance  curves . 

These  methods  provide  a  precise  procedure  for  criteria  development 
by  computing  exceedance  curves  with  a  mean  and  tolerance  limits  for 
all  possible  exceedance  curves  in  the  given  category.  The  results 
appear  promising  in  providing  a  more  descriptive  USAF  military  speci¬ 
fication,  which  in  turn,  necessarily  results  in  more  refined  structural 
design  for  future  aircraft. 
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SECTION  II 


DISCUSSION 

To  achieve;  the  objectives  stated  above,  several  methods  of 

analyzing  the  flight  loads  data  as  extracted  from  References  1  and 

2  were  employed.  For  each  set  of  data  by  aircraft  tail  numbers,  the 

following  statistical  functions  were  computed  for  the  vertical  load 

factor  (n  )  exceedance  curves: 
z 

a.  Maximum,  minimum  and  median  exceedance  curves  with  a  90% 
confidence  interval  for  the  median. 

b.  90%  two-sided  tolerance  limits  with  90%  confidence  assuming 
the  exceedances  are  normally  distributed. 

c.  95%  two-sided  tolerance  limits  with  95%  confidence  assuming 
the  log  of  the  exceedances  are  normally  distributed. 

d.  Least  squares  curve  fits  to  the  log  of  the  composite  data. 
Plots  of  the  above  results  are  shown  in  Figures  1  through  36  in 
Appendix  I  of  this  report.  Detailed  discussions  of  the  flight  data, 
computer  programs,  and  assumptions  used  to  obtain  the  above  results 
are  presented  in  the  following  paragraphs  of  this  report. 

1 .  Flight  Data 

The  load  factor  (n  )  peak  data  used  for  this  study  were  obtained 
z 

from  References  1  and  2.  Reference  1  contained  operational  flight 
loads  data  from  the  A-37B  aircraft  while  performing  training  flights 
from  May  1969  to  September  1971.  A  total  of  4,001  hours  of  data  was 
analyzed  and  is  summarized  in  Table  I.  For  a  detailed  description 
of  the  data,  the  reader  is  referred  to  Reference  l.t’ 


TABLE  I 


SUMMARY  OF  A-37B  DATA  USED^ 


Base  and  Year 

Hours 

Tail  Numbers 

Hours/Tail  Number 

England  AFB  (1969) 

541.47 

11 

49.22 

Bien  Hoa  AB  (1970) 

2038.23 

12 

169.85 

Binh  Thuy  AB  (1971) 

913.84 

8 

114.23 

England  AFB  (1971) 

507.87 

4,001.41 

4 

126.97 

TABLE  II 

SUMMARY  OF  F-4  DATA  USED^ 

Mission  Type  Hours  Tall  Numbers  Hours/Tall  Number 


Air-Ground 

2308.6 

46 

50.19 

Air-Air 

149.0 

23 

6.48 

In-Nav 

481.2 

40 

12.03 

Recon 

515.2 

19 

27.12 

Test 

19.1 

3,473.1 

15 

1.27 
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Reference  2  contained  F-4  aircraft  Southeast  Asia  (SEA)  load 
factor  data  that  were  segregated  by  tail  number,  and  mission  types. 
These  data  were  collected  by  Technology,  Incorporated  (TI),  as  part 
of  Aircraft  Structural  Integrity  Program  (ASIP)  during  the  period 
of  15  August  1969  through  31  December  1970.  The  3,473  total  hours 
analyzed  are  summarized  in  Table  II. 

The  above  referenced  reports  were  selected  since  the  data  are 
in  terms  of  aircraft  tall  number  by  mission  type  or  base.  Hov/ever, 
these  data  are  not  in  terms  of  mission  segment  as  desired.  It 
seems  reasonable  to  assume  that  segregation  by  base,  and  mission 
type  would  be  as  applicable  to  the  analysis  methods  as  segregation 
by  mission  segment. 

Before  proceeding,  some  of  the  problems  encountered  in  using 
the  above  data  should  be  considered.  One  problem  that  appeared  in 
Tables  I  and  II  was  that  both  aircraft  types  have  what  appeared  to 
be  a  large  total  number  of  hours,  but  may  have  too  small  a  number 
of  hours  for  any  given  tail  number  and  category  of  mission  or  base. 
For  example,  there  is  4,001  hours  of  F-4  data  but  only  an  average 
of  50  hours  per  tall  number  in  the  air-ground  mission  category,  with 
other  mission  categories  having  even  less  hours.  This  raised  the 
question  about  the  validity  of  an  exceedance  curve  for  a  given  tail 
number.  A  small  number  of  hours  for  a  tall  number  would  probably 
give  the  exceedance  curves  a  wider  spread  than  is  actually  the  case; 
this  would  tend  to  make  any  estimates  about  the  variability  to  be 
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unnecessarily  large.  If  the  criteria,  that  1,000  hours  of  data  are 
needed  to  plot  a  valid  exceedance  curve  is  used ,  it  would  be  necessary 
to  have  at  least  46(1,000)  +40(1,000)  +  19(1,000)  +  15(1,000)  — 

143,000  total  hours  of  F-4  data  to  obtain  valid  exceedance  curves 
for  all  the  aircraft  tall  numbers  and  categories  listed  in  Table  II. 
This  is  a  large  data  requirement. 

Another  question  concerns  the  type  of  distribution  the  exceedances 

have  for  a  given  value  of  n  .  It  will  be  shown  that  normality  should 

z 

probably  be  ruled  out,  but  the  log  of  the  exceedances  being  normally 
distributed  is  a  plausible  assumption.  The  assumption  of  a  nonpara- 
metric  (not  normal)  distribution  would  also  provide  some  usable 
results,  but  more  tail  numbers  would  be  required,  and  hence  more  data. 

Lastly,  the  problem  of  zero  exceedances  at  high  n^  values  caused 
a  number  of  interrelated  problems,  for  example  zero  peaks  in  50  hours 
implies  zero  peaks  per  1,000  hours  is  not  a  valid  assumption.  That 
is,  the  real  number  of  peaks  per  1,000  hours  at  some  high  load  factor 
may  have  really  been  10  peaks  implying  .5  peaks  in  50  hours.  Since 
there  can  be  no  fractional  count  of  a  peak,  0.5  peaks  in  50  hours  of 
flight  cannot  occur.  Even  with  over  1,000  hours  of  data,  the  result 
is  questionable  since  occurrences  of  n^  peaks  at  high  n^  values 
(above  7  g's)  are  rare.  Handling  of  this  zero  occurrences  problem 
is  discussed  throughout  this  report . 

It  was  decided  that  these  data  would  suffice  for  this  study, 
since  there  did  exist  variability  of  the  exceedance  curves  for  each 
aircraft  tail  number.  However,  for  better  results,  more  data  hours 
per  aircraft  and  more  tail  numbers  should  be  used. 


2.  Computer  Programs 


Several  computer  programs  (See  Figures  37-40)  were  used  on  the 

Hewlett-Packard  (HP)  9830  calculator  system  to  verify  the  validity 

of  some  new  approaches  to  analyzing  flight  loads  data.  The  computer 

programs  were  used  to  generate  plots,  tolerance  limits,  confidence 

intervals  and  curve  fits  for  the  n  exceedance  data  used  herein. 

z 

One  computer  program,  which  handles  a  maximum  of  49  tail  numbers 

and  16  different  values  of  n  along  the  abscissa,  used  the  n  peak  data 

z  z 

obtained  from  the  references,  and  stored  it  on  a  file  of  a  tape 

cassette.  Hence,  each  file  of  the  tape  cassette  contained  n  data 

z 

by  tall  number  for  a  particular  mission  or  base  as  shown  in  Tables 
I  and  II,  and  also  the  composite  data  which  includes  all  tall  numbers. 

The  tape  cassette  was  used  to  input  the  n^  data  in  ah  analysis 
program.  The  analysis  program  was  used  to  analyze  the  data  as 
follows : 

a.  Plot  n^  exceedance  points  per  1,000  hours  for  each  tail  number. 

b.  Plot  n  exceedance  curves  for  each  aircraft  tail  number. 

z 

c.  Compute  and  plot  the  minimum,  maximum  and  median  exceedance 

per  1,000  hours  for  each  of  the  16  n  values. 

z 

d.  Compute  and  plot  a  90%  confidence  interval  for  the  median  of 
the  exceedances  per  1,000  hours. 

e.  Compute  and  plot  the  mean  and  composite  exceedances  for  each 

of  the  16  n  values, 
z 

f.  Compute  and  plot  tolerance  limits  on  the  exceedances  per  1,000 
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hours  for  each  n  stratum  assuming  that  the  exceedances  or  the  log 
z 

of  the  exceedances  are  normally  distributed. 

A  third  computer  program,  which  was  a  revision  of  a  HP  Plot  Pac 
Program,  performs  a  least  square. fit  of  a  polynomial  to  the  log  of  the 
exceedance  data.  Composite  data  for  all  tail  numbers  for  each  category 
were  input  through  the  keyboard  to  fit  curves  of  the  form 

F  . 

where  F  is  the  number  of  exceedances  per  4,000  hours  and  b(n  )  is 

z 

the  best  polynomial  that  estimates  log  F. 

3.  Data  Presentation 

Results  of  the  analysis  are  presented  in  Figures  1  through  36  and 

Table  III  in  Appendix  I.  A  detailed  discussion  of  the  various  methods 

! 

that  were  studied  is  presented  below. 

3.1  Exceedance  Curves  by  Aircraft  Tail  Numbers 

For  each  of  the  nine  categories  listed  in  Tables  I  and  II, 
the  exceedance  curve  for  each  individual  tail  number  has  been  plotted 
in  Figures  1  through  9.  These  exceedance  curves  give  a  rough  idea 
of  the  type  of  spread  that  might  be  expected  for  different  aircraft 
tail  numbers.  The  following  sections  primarily  investigate  the  means 
by  vdiich  statistical  statements  can  be  made  in  describing  this  spread 
in  the  exceedance  curves . 

3.2  Confidence  Intervals  for  the  Median 

The  median  of  the  exceedances  per  1,000  hours  (MED(n  ))  for 

Z 
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each  given  value  of  have  been  plotted  in  Figures  10  through  18 

for  each  of  the  nine  cases  listed  in  Tables  I  and  II.  Straight 

line  segments  were  also  drawn  connecting  the  MED(n^)  values.  The 

median  was  calculated  so  that  one  can  say  that  about  1/2  of  the 

exceedances  per  1,000  hours  for  each  n  (E.  (n  )  where  1  is  the  tail 

z  1  z 

number)  lie  above  and  below  the  MED(n  )  value. 

z 

Next,  a  90%  two-sided  confidence  interval  (Cl)  for  the 

median  was  found  using  a  distribution  free  procedure  based  on  the 

sign  test  as  outlined  in  Reference  3.  The  procedure  was  applied  to 

each  E.(n  )  for  a  given  n  so  that  the  90%  confidence  interval  for 
1  z  z 

the  median  corresponds  to  the  statement:  the  probability 

LCI(n  )  <_  MED(n  )  <.  UCI(n  )  =  0.90 

where  LCI(n^)  is  the  lower  limit  of  the  Cl  and  UCI(n^)  is  the  upper 

limit  of  the  Cl  (i.e.  the  probability  that  the  true  median  lies 

between  LCI(n  )  and  UCI(n  )  is  0.90). 
z  z 

All  E^(n^)  =  0  cases  were  included  for  the  above  computation. 
In  addition,  the  composite  was  also  plotted  and  found  within  the  90% 

Cl  for  the  median  for  the  lower  values  of  n^.  Figures  10  through  18 
showed  that  a  difference  exists  between  the  median  and  the  composite 
(C)  for  the  F-4  data  whereas  they  are  very  similar  for  the  A-37  data. 
In  fact,  for  the  F-4  data  the  median  usually  lies  below  the  composite. 
This  is  probably  due  to  the  greater  frequency  of  E^(n^)  *  0,  for  the 
F-4  data  which  was  caused  by  an  inadequate  number  of  hours  per  tail 
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number  as  shown  in  Table  II  and  discussed  earlier  in  data  presenta¬ 
tion.  It  should  be  noted  that  1  peak  in  100  hours  of  data  would 
imply  10  peaks  pet  thousand  and  may  overestimate  a  theoretical  value 
of  1  peak  per  thousand  hours  by  quite  a  margin.  Intuitively ,  it  then 
seems  that  for  the  F-4,  the  composite  would  usually  lie  below  any 
exceedaiice  per  tail  number  due  to  the  overestimation  of  the  exceedance 
for  a  limited  number  of  hours. 

3.3  Maximtan  and  Minimum  Values 

The  maximum  and  minimum  values"  of  Ej|^(n^)  are  also  plotted  in 

Figures  10  through  18.  The  maximum  (MAX(n^))  and  minimum  values 

(MIN(n  ))  are  plotted  for  each  n  value  and  a  line  drawn  connecting 
z  z 

the  points.  These  curves  may  then  be  used  to  make  distribution  free 
two-sided  tolerance  limit  statements  which  means  that  with  some  confi¬ 
dence  gamma  (y) «  a  percentage  of  the  tall  numbers  (P)  lies  between  the 
two  limit  values.  By  using  Table  A-32  in  Reference  4,  values  of  y 
are  obtained  so  that  90%  of  the  tail  numbers  lie  between  the  maximum 

and  minimum  values  of  n  exceedances  per  1,000  hours. 

z 

In  Figures  10  through  18,  the  values  of  y  have  been  included 
on  the  plots.  Note  that  gamma  was  higher  as  the  number  of  tail  numbers 
increased.  For  example,  y  =  .30  for  11  tail  numbers  in  Figure  10 
and  y  =  0.95  for  46  tail  numbers  in  Figure  14.  Hence,  we  can  use  the 
maximum  and  minimum  value  of  exceedances  per  1,000  hours  as  statistical 
bounds.  This  method  would  give  a  very  conservative  estimate  of  the 
bounds  on  the  exceedance  curves  since  there  are  no  assumptions  made 


about  the  underlying  distribution.  For  Figures  10  through  14,  the 
bounds  determined  in  this  manner  appear  to  be  relatively  narrow. 
However,  in  Figures  15  through  18,  the  minimum  is  zero  too  often 
since  these  cases  have  a  small  average  hours  per  tail  number. 

3.4  Tolerance  Limits  Assuming  Normal  Distribution 

In  Figures  19  through  27,  a  plot  of  a  90%  two-sided  toler¬ 
ance  limit  with  90%  confidence  is  shown  with  the  assumption  that  for 
a  given  n^,  the  exceedances  per  1,000  hours  have  a  normal  distribu¬ 
tion. 

The  tolerance  limit  was  computed  by  considering  the  E  (n  )s 

i  z 

for  each  given  n^.  The  mean  M(n^)  was  calculated  using 


M(n  )  =  2  E,  (n) 

^  i=l  -  , 

N 

where  (n  )  is  the  exceedances  per  1,000  hours  for  a  fixed  n  and  N 
-*-2  2 

is  the  number  of  aircraft  tall  numbers.  Then,  the  sample  standard 
deviation(s)  were  computed  by 


s(n^)  = 


n  E  E  (n  -  (  Z  E  (n 

i"l  ^  ^  i=l  ^  M 


n(n-l) 


which  was  obtained  from  Reference  1  on  page  1-10. 
The  upper  NorUL  and  lower  limits  NorLL  were 


NorUL(n  )  =  M(n  )  +  K*s  (n  ) 
z  z  z 
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NorLL  =  M(n^)  -  K‘S  (n^) 

where  K  is  obtained  from  Table  A-6  of  Reference  4  with  y  =  .90, 

P  “  90%  and  N  “  number  of  tail  numbers.  Before  proceeding,  it  should 
be  noted  that  all  zero  exceedances  per  1,000  hours  were  included  in 
the  computations . 

The  lower  limit  did  not  appear  in  most  of  the  figures , 
because  the  assumption,  that  the  exceedances  are  normally  distributed, 
allows  the  lower  limit  to  take  on  negative  values.  This  is  not 
possible.  Another  factor  to  consider  concerns  plotting  of  limits 
on  a  log  scale.  Since 

1  log  (M  +  KS)  -  log  M  I  <  1  log  (M-KS)-log  M  | 

the  upper  limits  would  be  closer  to  the  mean  than  the  lower  limits 
on  the  log  scale. 

As  was  the  case  for  the  median,  the  mean  and  composite  are 

approximately  equal  at  low  n^  values  but  differ  at  higher  n^  values. 

This  is  again  probably  due  to  errors  in  exceedances  per  1,000  hours 

that  occur  at  the  higher  n  values. 

z 

It  was  concluded  from  the  above  discussions  and  Figures  19 
through  27  that  the  assumption  of  a  normal  distribution  was  not  valid. 
3.5  Tolerance  Limits  Assuming  a  Log  Normal  Distribution 

Plots  of  90%  two-sided  tolerance  limits  with  90%  confidence 
are  shown  in  Figures  19  through  27  with  the  assumptions  that  the  log 


of  the  exceedances  per  1,000  hours  for  each  has  a  normal  distri¬ 
bution.  The  mean  M3  was  computed  as  the  common  logarithm  of  the 

composite  (C(n  ))  at  each  n  .  Note,  the  composite  is  the  same  as 
z  z 

the  weighted  mean  (M^)  defined  by 


V  E.(n  ) 

M  =  E  .1  z 

”  i=l  Et^ 

which  equals  the  total  number  of  occurrences  greater  than  n  divided 

z 

by  total  time.  This  result  is  the  composite.  Next,  a  standard 


deviation  s'(n  )  is  defined  as 
z 


Since,  the  log  of  zero  is  undefined,  zero  values  of  E . (n  )  were  not 

1  z 

included  in  the  solution  for  s'(n  ).  N  was  decreased  by  1  each  time 

z 

a  zero  was  excluded  for  a  fixed  value  of  n  . 

z 

Then,  as  in  the  previous  section,  upper  (LNUL)  and  lov/er  (LNLL) 
tolerance  limits  were  computed  by  formulas: 


LNUL(n,)  =  log  C(n^)  =  K-S'(n^) 
z  z  z 

and 

LNLL(n  )  =  log  C(n  )  -  K*s'(n  ). 

Z  2  Z 

These  upper  and  lower  tolerance  limit  (LNUL,  LNLL)  curves  are  then 
equidistant  from  the  composite  when  plotted  on  semilog  paper  as  is 
evident  from  Figures  19  through  27.  For  most  cases,  the  log  normal 
assumption  gave  some  reasonable  bounds  but  deviations  do  occur.  For 
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instance,  at  the  extreme  n  values,  the  limits  go  to  infinity  and 

.  z 

zero  for  the  upper  and  lower  bounds  respectively.  This  is  because 

at  the  extreme  n  values  there  are  many  zeroes  which  make  N  small , 

■  .z  '  ' 

and  K  which  was  obtained  from  Table  A-6  becomes  very  large.  Hence, 
LNUL  and  LNLL  become  very  large  and  small  respectively  as  can  be  seen 
clearly  in  Figure  26.  The  spike  at  n^  -  6  g’s  on  Figure  22  resulted 
from  one  aircraft  tail  number  exceedance  curve  being  farther  from 
the  mean  than  the  other  three;  thus,  making  a  higher  than  probable 
standard  deviation.  It  was  decided,  that  overall,  the  log  normal 
assumption  did  provide  reasonable  bounds  for  n^  exceedance  curves. 

3.6  Composite  Exceedance  Curve  Fitting 

The  final  computer  program  provided  for  the  fitting  of 
polynomial  equations  to  the  log  of  the  exceedances  per  4,000  hours. 
Four  thousand  hours  was  choseli  because  it  is  the  design  life  for  most 
fighter-type  aircraft.  The  fitted  curves  that  were  obtained  are 
presented  in  Figures  28  through  36  and  summarized  in  Table  III. 

In  most  instances,  a  third  degree  polynomial  provided  the 
best  fit,  but  in  instances  where  this  was  not  a  good  fit,  (e.g.  the 
curve  goes  to  plus  infinity  as  n  goes  to  infinity)  the  best  straight 
line  was  fitted  through  the  data  points.  In  most  cases  when  the 
straight  line  fit  was  used  it  proved  to  be  adequate. 

These  fitted  curves  to  the  n^  exceedances  per  4,000  liours 
can  be  used  as  a  uniform  method  of  extrapolation  to  one  exceedance 
per  4,000  hours  (i.e.  the  aircraft  life  to  obtain  an  estimate  of  a 
realistic  design  limit  load  factor). 
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In  addition,  these  equations  may  be  used  to  change  the 
magnitude  and  slope  of  the  exceedance  curves  for  test  purposes . 

The  equation  has  the  form 

F(x)  = 

where  F  is  the  curve  function,  x  corresponds  to  n  and  B  is  a  function 

z 

of  X  of  the  form 

B(x)  =  +  bj^x  +  Bq 

For  the  case  of  the  cubic  equation  for  positive  n^  in  Figure  32,  the 
equation  for  the  exceedance  curve  for  F-4  air-ground  data  was  found 
to  be 

(6. 9742-1. 4829x  +  . 2914x^-0. 0249x^) 

F(x)  -  10 

(See  Table  II  in  Appendix  I) 

This  equation  can  be  reduced  to 

, -6.9742, _x(-.0249x  +  .2914x-  0.0249x^) 

F(x)  =  10  10 

This  means  we  can  place 

P  =  10^ •9742^  9,423,234 

and 

B(x)  -  (-.0249x^  +  ,2914x-1.4829) 

P  will  define  the  point  at  which  the  curve  intercepts  the  n^  =  0 
axis  and  B(x)  defines  an  instantaneous  slope  to  the  exceedance  curve 
on  the  log  scale.  This  takes  into  account  the  changes  in  slope  which 


is  typical  for  fighter  exceedance  curves.  As  long  as  is  a 
negative  coefficient  (-.0249  for  the  example  cited  above),  we  know 
that  the  fitted  curve  will  approach  zero  as  n^  Increases  without 

bounds. 

With  these  equations,  a  series  of  curves  could  be  obtained 
by  varying  the  magnitude  P  and  coefficients  of  the  slope  B(x)  such 
that  they  do  not  exceed  the  bounds  determined  by  using  one  of  the 
previous  techniques.  The  effect  of  varying  the  P  and  the  B(x)  has 
not  been  studied  in  this  report,  but  might  be  an  approach  that  may 
prove  effective. 

Extrapolation  of  these  curves  to  design  limit  load  factor 
is  now  a  simple  procedure.  Since  F  is  the  exceedances  per  4,000 
hours,  find  the  value  of  n^  such  that  the  log  F  =  log  1  =  0.  For 
example  in  Figure  28  for  A-37B,  England  AFB,  1969,  we  would  find 
a  design  limit  load  factor  of  6.8  g's. 
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SECTION  III 
CONCLUSIONS 


The  plotting  of  exceedance  curves  by  aircraft  tail  number 
to  obtain  statistical  bounds  did  appear  to  be  feasible.  Finding 
two-sided  tolerance  limits  for  exceedance  curves  with  the  assumption 
that  the  log  of  the  exceedances  per  1,000  hours  are  normally  distri¬ 
buted  appeared  to  give  reasonable  bounds.  Nonparametric  techniques 
gave  reasonable,  but  more  conservative  bounds.  The  assumption  that 
the  exceedances  per  1,000  hours  were  normally  distributed  appeared  to 
give  the  widest  bound  and  hence  were  not  as  useful .  Curve  fitting 
polynomials  to  the  log  of  the  exceedances  per  4,000  hours  for  compo¬ 
site  n^  data  provided  usable  equations  for  representing  n^  exceedance 
curves  in  most  cases  for  the  two  aircraft  types  that  were  considered. 
In  addition,  the  polynomials  provide  a  standard  mean  of  extrapolating 
n^  exceedance  curves  to  design  limit  load  factor. 


SECTION  IV 


RECOMMENDATIONS 

As  noted  earlier,  the  data  used  for  the  F-4  represented  a  small 
amount  of  flight  time  per  tail  number.  A  further  study  should  be 
made  to  determine  what  minimum  time  on  a  tail  number  is  required  in 
each  category  to  ensure  that  a  representative  exceedance  curve  can 
be  obtained  for  that  particular  aircraft.  The  A-37B  data  indicated 
that  100  hours  per  aircraft  tail  number  would  probably  be  sufficient, 
but  more  Information  would  have  to  be  studied  before  this  statement 
could  be  made  with  certainty. 

Separating  the  data  into  mission  segments  would  Increase  the 
need  for  a  larger  data  base.  For  example,  if  nonparametric  techniques 
are  used,  then  at  least  50  tail  numbers  should  be  used  for  each 
mission  segment.  If  100  hours  per  tail  number  are  required  and  five 
mission  segments  for  a  typical  aircraft  type  are  assumed,  there  would 
need  to  be  at  least  25,000  hours  of  data  on  the  particular  aircraft 
type.  If  the  data  were  not  categorized  by  mission  segment,  5,000 
hours  would  suffice,  and  hence,  it  is  seen  that  the  total  number  of 
hours  needed  is  directly  proportional  to  the  number  of  categories 
within  the  data  base.  Thus,  it  should  be  remembered  that  segregating 
data  into  mission  segment  would  greatly  increase  the  amount  of  data 
that  needs  to  be  collected. 

Also,  the  uncertainty  of  the  nvimber  of  exceedances  per  1,000  hours 
at  higher  n^  values  is  of  concern.  A  minimum  of  100  hours  per  tail 
number  would  be  a  minimum  requirement;  however ,  there  is  still  a  problem 
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when  zero  exceedances  per  1,000  hours  occur,  because  exceedance  per 
1,000  hours  approaches  zero  as  increases,  but  does  not  really 
reach  zero.  These  problems  at  high  n^  values  could  easily  form  the 
major  topic  of  further  study. 

Lastly,  the  aircraft  for  this  report  were  both  fighter-type 
aircraft.  For  other  aircraft,  like  transport,  other  techniques  may 
be  more  applicable  than  those  cited  in  this  study. 
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APPENDIX  I 
DATA  RESULTS 
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5.  Exceedance  Curves  by  Tail  Number  for  F-4  Air-Ground 
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FIGURE  6.  N  Exceedance  Curves  by  Tail  Number  for  F-A  Air-Air 
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FIGURE  10.  Nonparametrlc  Bounds  for  A-37B  at  England  AFB 
(1969)  Exceedances 
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FIGURE  11.  Nonparametric  Bounds  for  A-37B  at  Bien  Hoa  AB 
(1970)  Exceedances 
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Confidence  that  90%  of  tail 
numbers  lie  between  the  maximum 
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FIGURE  12.  Nonparamet’rlc  Bounds  for  A-37B  at  Binh  Thuy  AB 
(1971)  Exceedances 
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FIGURE  13.  Nonparametrlc  Bounds  for  A-37B  at  England  AFB 
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FIGURE  15.  Nonparametrlc  Bounds  for  F-4  Air-Air  Exceedances 
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FIGURE  17 .  Nonparametric  Bounds  for  F-4  Reconnaissance  Exceedances 
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FIGURE  18.  Nonparametric  Bounds  for  F-4  Test  Exceedances 
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FIGURE  19.  90%  Tolerance  Limits  with  90%  Confidence  for  A-37B  at 

England  AFB  (1969) . 
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FIGURE  21.  90%  Tolerance  Limits  with  90%  Confidence  for  A-37B 

at  Binh  Thuy  AB  (1971) . 
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FIGURE  22.  90%  Tolerance  Limits  with  90%  Confidence  for  A-37B 

at  England  AFB  (1971). 
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-2  -1  0  1  2  3  4  5  6  T' 

FIGURE  23.  90%  Tolerance  Limits  with  90%  Confidence  for  F-4 

Air-Ground. 


8  n^(g' 
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40  Tail  Numbers 
481.2  Total  Hours 
ASD-TR-71-37 


I  15  Tail  Numbers 
j  19.1  Total  Hours 
'  ASD-TR-71-37 


I  i  I  I  1  -  Upper  bound  of  tolerance  limit  assuming  normal 

!  ;  ;  ;  2  -  Mean 

i  ;  —  ^  ~  Lower  bound  of  tolerance  limit  assuming  normal 

ri  I  .  4  -  Upper  bound  of  tolerance  limit  assuming  lognarms 

;  1  i  ^  5  -  Composite 

I  I  I  ;  6  -  Lower  bound  of  tolerance  limit  assuming  lognarms 


-1  0 


8  n^(g's) 


FIGURE  27.  90%  Tolerance  Limits  with  90%  Confidence  for  F-4  Test. 
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TABLE  III.  CURVE  FIT  COEFFICIENTS 
3  2 

F(X)=10^®3^  “^^2^  ■'■®l^'*'®o^where  x=n  and  F(X)=  Exceedances /4 ,000  Hours 

z 

A.  Coefficients  for  n^  >lg 


Figure 

!3 

fl 

a 

0 

R-Square 

Hours 

28 

-0.0592 

.5281 

-1.9186 

7.3527 

.99908 

541.47 

29 

-0.0305 

.2907 

-1.3611 

6.6859 

.99796 

2038.23 

30 

-0.0954 

.8393 

-2.7681 

7.6166 

.99956 

913.84 

31 

-0.0254 

.2390 

-1.1635 

6.9181 

.99913 

507.87 

32 

-0.0249 

.2914 

-1.4829 

6.9742 

.99681 

2308.6 

33 

-0.0177 

.1833 

-1.0033 

6.4857 

.99945 

149.0 

34 

-0.0525 

.5726 

-2.6019 

7.7763 

.99583 

401.2 

35 

-.6364 

6.4918 

.98916 

515.2 

36 

-0.0117 

.1180 

-.8723 

6.3584 

.99442 

19.1 

B.  Coefficients 

for  n  < Ig 
z 

Figure 

!3 

!2 

!l 

a 

0 

R-Square 

Hours 

28 

3.4216 

4.2709 

.99962 

541.47 

29 

.6777 

2.4887 

4.4009 

4.1837 

1 

2038.23 

30 

1.7549 

3.3475 

.87189 

913.84 

31 

2.4042 

4.4226 

.99650 

507.87 

32 

.2336 

.9398 

2.4319 

2.9379 

.99968 

2308.6 

33 

1.6055 

3.7242 

.98488 

149.0 

34 

2.7188 

2.7525 

.99223 

401.2 

35 

2.5057 

3.2584 

.99435 

515.2 

36 

.5177 

1.2371 

1.2208 

3.5295 

.99835 

19.1 
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FIGURE  29.  Composite  Data 


N2(5} 

FIGURE  30.  Composite  Data  Curve  Fit  for  A~37B  at  Binh  Thuy  AB  (1971) 


/  OS) 

FIGURE  31.  Composite  Data  Curve  Fit  for  A-37B  at  England  AFB  (1971) 


FIGURE  32,  Composd 


'-4  Inst,  and  Nav. 


APPENDIX  II 


COMPUTER  PROGRAMS  USED 


FIGURE  37.  Program  to  Store  Data  on  Tape  Cassette 


;  iFiJi-!  Fi0 

[  50  ]:■  NSil  50?  16  ]  ? 

■  r  F;R  J- 

1  T'O  150 

::3  H[J1=0 

1  TO  16 

A  :  I  1  '■ 

ip  t[J]=^m;j 
:p  vlp. 

....  {.J 

.7  HE  NT  l< 

.10  IPiPp  " 

R0WL^==?i:;i:iLUi'1NM=" 

iiiPU! 

L..  ?  I'1 

■'0  FOR  I=- 

1  TO  1... 

I'O  IJISP  “ 

FK  '2[  '  ):=?" ; 

;10  INPUT 

pi:  I  ] 

7  0  FOR  ..P= 

1  TO  M 

70  DISP  " 

N(  “  I  ?  ■'  ?  "  J"  )=?''  ii 

00  INPUT 

N[  I  ?  J  ] 

NENT 

.J 

110  BISP 

"  T  C "  I "  ? 

,120  INPUT 

Ti;  I  ] 

130  NENT 

I 

170  FOR  J 

=  i  TO  11 

156  DISP 

”  Y  (  "  J "  )  =  ”  5 

leO  INPUT 

YC  ..,,1  ] 

170  j 

1S0  FOR  J=1  TO  L 

100  PRINT 

;7  00  prihT 

:l;i6  FOR  K==l  TO  8 


VSC  16] 


PRINT 

FOR  TO  8 
100 1  NT  V[  I  ]; 
HE  NT  I 
PRIHT 


277  PRINT 

200  FOR  J=1  TO  L 

205  FOR  10" 9  .  TO  16 


FI6U§^  37.  Continued 


Ui-iV  XL  J  5  K  j  :■ 

:.;:';r  k 

'••:.:  J 

■•■::  '■'L-Nl' 

:'iJK  ’.-■y  MJ  16 

l■■^■.iHT  VLIj; 

■'4S  I 

5  FRIHT 

I' ISP  ''SNTER  1  :F  1X17 

i:r:riJ7  ti 

'  1  i  ^  THEN  4  40  . 

Sv--;!  ‘DHIH  Bfil.i  i'H  !  fi 

xyy  SPLIT 

‘  j  L-  0_  .!  .i  y  F  '  i -i  L  "Li?  "  v;  ( '■  j  :  ■■ 

4 1  p:  N P I, i T  F! U..:  1  ] .  Li L  L.,  I  ,  M 1  ii  j  'f  li  L  U  iF/r  M 1  ] 

‘  ‘  !'■  ill  i  i  :[.  L  I  J  :i  r,!.  L.  J.  J  f1  ,i.  J  5  i  L  L  1  j  I*  Y!  i'l  .1  J 

4:40  SOTO  ,;.:;5o 

4  -1  41  SI  '  ”riLi:'  NO„r." ; 

S'  '4  i'ilL-OT  z 

■■'"'■l  :4T0R,':;  IlFlTR  Z 
4  Tup 
4L:0  end 


,  0  bSiOD"  ? 


j  ?  T  ( “  L 1 5  “  >  Y  (  '■ 


FIGURE  38.  Program  to  Read  Data  on  Tape  Cassette 
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FIGURE  39 .  Continued 


r-, 


•i  I 


'Jr  -i  -rp 

■Pirr  Ti 


'p  JJJfl  then  in 00 
ui< 


.,.1 


i'  l:N  !■' ! 


■  0  r  iJJNT  Ri 


'yr  -.N  SH  "ENTER  1  ro  PR IRT  '  Cr' S "  5  ■ 
■ir  J:-iPUi  Ti 

J:-;::''  1390 

!'■■!  L.  I-j;  ;  \i  i 
7’  I:;' 

■■ '  '  j' : ! ' !  i .. ).  I .  U 1'.!  i:.  f'.;  J  Q  [ .  [  _  1 j;  j  |  p  p  r-  r;. 

rj-  '  J . enter  i,  TO  FlNJi  ri?,'vp%'cp 

JnENT  Ti  ■  . . 

Jr  N'  i  1  Trj,  THEil  1816 
:■  TE.R  TO  N 
r.'  J'-J'N'  i  ■:  IJ  J 

T  'J  ■  '=r'  T  Jp  K  Ti'Q 


...  ...| .  ^ 


i'OR  .i  =  S^  Ti'i 


N 

IMG  NORitni 


t- 

I 


f 

[ 

f 


( 


Ft§URE  39.  Continued 


I  .  i  i.'i  ! 


’010  nSF'  "ENTER  i  TO  FIND  ME ■ ’^^0 ?  S2 ^  HSS- 'h  f  HE  lOE  NORMR 


l!  I  ST  R I BU  ■ 


■■E0  lilPuT  Ti 


'1  j  r  LI 


040  '  I'. T-:  K>"1  TO  i1 
O'::'0  IF  Fi: ...  ? J-TT  i'MiiH  -El  EO 
i "-Lij  i  HI  I...  >  i  -.  j 
;  .j™  ,!.  1  U  N  , 

,  ■’30  tF'  xi:  ti-en  jf/Eo 

.  ■?  0  0  “  ( L  .  G  T  X  i :  J  5  R:  j  -  P'i  E  )  ■  £ 

1919  GOTO  193.0 

:;.9£0  i!EX'!'  J 

■;  ■■*30  FOR  ...l:-J-i-l  TO  H 

’iSTO  Jr  F;[J?E]=0  THEN  1980 

.!.  4  5  01  0 1  ~  i ...  Li  T !  i  [  -..i  ii  K  1  -- 1'  1  £  ) '  F  2 + U 

1 9£0  N  =  -L’i  '  ■ 

::  !.:.i  ..3  1  4  9  O 

998  T^iKX'i  J 

:::880  ir  n:9==i  then  Ftee 

£':il0  .1  F  ri'l!:;.  !■!  FlbM  '■ 

.i'F'SP  ‘  ER'-TfER  TO  Fll-if'  l  i  f  !..£■■  LOG  HOFHNL  BISTR,,  '5 


T£T  1  THEi-i  .£  i.  TO 


Mi- ...I  i 


F  i  00  L  .9=4'E24K5-v3£ 

L 1 1 H  r  r  1 ,  !<  ]="■•/ 2 


3130  C:-3rE]:442 

iL!  C*  rO  !  '  *  !:L  E  l  I  R  . 

103  :;NPi.iT  T.i. 

1  !'■  Tiiti  THE;^^; 

.  i;3.  3 1--,.':  T*)  -■ 


!  U  I’iJJ  i  L  i  f  Me  f  LE  > 


i  1 .  i  3  3  ‘  ..i  i '' 
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FIGHSe  39 .  Continued 


-!PiJ 


30 


i;;} 


t- Pi 00  thpn 

0'  i'i  ]  :■  iji"  13  1  :i  [:;■  ] 
i' THEN  ;;:':4'30 


Pi.0  i'EM 
i  0  EE  r;Jj3N 


Ei''  '.M-ii  Ci"  K  1  J  J  J 

0  iiETijRH 
is  iE'-ill 
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# 


J 

FIGURE  40.  Continued 


10  IF  W  THEIN  60 

15  niEiP  "KHOilpFFj 

16  IF PUT  T 

20  DISP  "XjV" 


30  IIIPIJT  BI;21hY 
32  PfllHT  "HZ-"BC2 
o4  V=V/T 


CUMULRTIVE  PERKS 


:J5  'T  =  Y:':'40Ri"i 
36  PRINT  "NZ 
'  I'  r  L.  l:i  1  Y 


BE 


CUMULRTIVE  PERKS 


'10  IF  FNNl  THEN  SO 
30  END 


60  DISP  "NOT  RLLOWED 
70  END 


PER  4000  HR 


0  IF  N  THEN  60 


20  DISP  "WRONR  FLY ^3' 
30  INPUT  B[2]?Y 
40  IF  FNN'::-!)  THEN  2 
50  END 


|l 


0 


68  :DIE::P  "NOT  RLLOHED" 
?'0  END 


■  1 E'  FI  9  =  ( S  5 

-  S 1  *  S  3  H )  <  N  - 1 )  /  S  8  /  S  9 

40  print" 

50  PRINT 

'NO.  POINTS  =“N 

^00  PRINT 

70  PRINT 

A  •  I'll;;. H N  =  '  S 1 ..' N 1  |  H B 2 .5  ' 

ST,DEV.=" SR 

00  PRINT 

''t'!  NEflN=“S3--'N;  TRB25‘ 

ST.DEV.="S9 

90  PRINT 

100  PRINT 

"CORR'.  COEFF.  =  "R9 

110  PRINT 

; 

120  END 

10  IF  N  <-  D2-W  THEN  260 
2  0  DISP  "D£G.REia.="; 

30  INPUT  D1 

40  IF  D1  <=  DS-N  THEN  70 
bu  DISP  "MRN  IiEG=";d£-h 
60  END 

70  IF  H=0  THEN  250 
00  T  =  0 

90  FOR  1=1  TO  Dl+1 
100  BE  I  ]=0 

110  FOR  ...1=1  TO  D1 . 1  +  R 

1 2'  0  F:  =  <  I  +  ...I  - 1  > ( j]  £  2  @ .  5  .JK  I  ■+•  J)) 
130  BE  I  ]=BE  I  ;i+CCT+J3*r:EP3 
140  NENT  J 

■i  so-  T=T-f  (■TlP-l-CR-T  2 > 


66 


^  -  p 


FIGURE  40 .  Continued 


160  NEKT  I 
170  R1=0 

180  FOR  1=2  TO  Dl+1 
190  R1=R1+CC  r*(ri2+<3-I>/2)  ]t2 
200  NEXT  I 

2 1 0  T0=CC  C  D2+ 1 ) * ( D2+2 ) /2  3 
220  T0=T0“CCD2+nt2 
230  DISP  "DONE" 

240  END 

250  If  N>D2  THEN  280 
260  DISP  "NOT  ENOUGH  POINTS" 

270  END 
280  P=W=1 
290  D2=D2+1 
300  FOR  J=1  TO  D2 
310  CCP3=SQRCCP3 
320  FOR  1=1  TO  D2-J+1 
330  CCP+I]=C[P+n/CCP-] 

340-  NEXT  I 
350  R=P+I 
360  S=R 

370  FOR  L=1  TO  D2-J 
380  P=P+1 

390  FOR  N=1  TO. D2+2-J-L 
400  C:[  R+M- 1 3=i;:C  R+N~  1  3-C[  P  ]*Cr  P+M**  1 3 
410  NEXT  M- 
420  R=R+M-1 
430  NEXT  L 
440  P=S 
450  NEXT  J 

4  6  0  1  =  '•  D  2  + 1  *  (  D  2 + 2 )  /  2 
470  FOR  1=1  TO  D2-1 
480  T=T-1-I 
490  C[T]=1/C.CT3' 

500  FOR  J=1  TO  D2-I 
510  P=D2+1-I-J 
520  P=P*(D2+l-<p-l>/2>-I 
530  R=P-J  ■ 

540  S=0 
550  iJ  =  I+.J+l 
560  V=P 

570  fOR  K=1  TO  J 

580  V=V+U-K 

590  S=S~CC  R+K  3*CC  V  3 

600  NEXT  K 

610  CCP3=S/CCR] 

620  NEXT  J 
630  NEXT  I 
640  CC  1  3=1/CC  13 
650  GOTO  80 


V 


FIGURE  40.  Continued 


IF  |.’.|  =  0  THEN  120 
PR  I  NT 

I  RIHT  "COEFFICIEHTS" 

J  PRINT 

R  FOP  IT  A  T  F3.0iiF12.4 
N  FOR  1  =  1  TO  Iil-i-l 
' !  ^  IT  R 1 1"  E  (15  ?  5  0 ) "  B  < " '[  1 "  J  = "  B  f  I  3- 
.0  NENT  1 
■hT  PRINT 

00  PRINT  "R  SQIJhRE  =  "R1/T0 
.10  PRINT 
.20  END 


0  !■  UR  K=N1  TO  N2  STEP  (K2'-Kl)/i00 

' .  *  I  . I  I  ‘i  C-  A 

■Ri  IF  'T<Y5  OR  Y>ViS  THEN  SO 
0  PLOT  NjY 
0  GOTO  70 
0  FNLN 
0  NENT  N 
0  Z=FNi<0 
0  END 


0  D I  P  '■  C H 13 F:  fl  C  T  E  R  H  E I G H  T  f  Z ) 
0  INPUT  H 

!....  H I:::' E L  (,  #  :i  H  j  2  ’  U-’  2/ Pi  ) 

0  LETTER 
0  Z-FNK0 
0  END 


DISP 
0  INPU 
,i  nisp 
0  END 


C) 


X 
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